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INTRODUCTION 


Breast  cancer  is  the  most  frequently  diagnosed  cancer  and  the  second  leading  cause  of  can¬ 
cer  mortality  among  women  in  America  [1],  Researchers  are  actively  developing  breast  imaging 
techniques,  including  breast  tomosynthesis  [2]  and  dedicated  breast  computer  tomography  CT 
[3,4],  that  can  form  three-dimensional  (3D)  images  of  the  breast,  thus  potentially  overcoming 
inherent  limitations  of  conventional  two-dimensional  (2D)  mammography.  Several  prototypes  of 
dedicated  breast  CT  scanners  have  been  built  in  different  institutions  throughout  US.  Current 
scanners  employ  a  circular  imaging  configuration  for  data  collection  and  the  FDK  algorithm  for 
approximate  image  reconstruction.  However,  it  is  well-known  that  a  circular  imaging  configuration 
does  not  yield  sufficient  data  for  mathematically  exact  image  reconstruction.  Also,  because  the 
FDK  algorithm  requires  data  collected  at  a  large  number  of  projection  views,  and  because  the  to¬ 
tal  dose  in  breast  CT  imaging  is  typically  about  the  same  as  that  in  a  two-view  mammography,  the 
number  of  photons  collected  at  each  view  can  be  severely  limited.  As  such,  the  projection  data  at 
each  view  have  very  low  signal-to-noise  ratio  (SNR),  thus  severely  limiting  the  breast-CT-image 
quality. 

Clearly,  image  quality  and  radiation  dose  are  great  concerns  in  breast  CT  imaging.  The  ob¬ 
jective  of  this  research  is  to  investigate  and  develop  innovative  imaging  configurations  and  recon¬ 
struction  algorithms  for  obtaining  accurate  images  and  reducing  radiation  dose  in  breast  CT  imag¬ 
ing.  In  the  past  year,  my  efforts  for  the  project  have  been  supported  by  a  Predoctoral  Traineeship 
Award,  and  I  have  performed  research  tasks  as  planned.  As  discussed  below,  I  have  studied  and 
assessed  several  potentially  useful  non-circular  configurations  for  breast  CT  imaging  by  investi¬ 
gating  data  conditioning  techniques  and  reconstructing  images  for  the  configurations.  I  have  also 
conducted  preliminary  studies  using  innovative  optimization-based  algorithms  for  a  potentially 
substantial  reduction  of  radiation  dose  in  breast  CT  imaging.  In  addition  to  computer-simulation 
studies,  I  have  performed  real  data  studies  involving  physical  phantoms  and  patients  provided 
by  the  collaborators.  In  summary,  during  the  first  year,  I  have  successfully  carried  out  research 
on  the  tasks  planned,  and  the  results  obtained  have  formed  a  solid  basis  for  me  to  continue  the 
research  planned  for  the  next  year. 


BODY 


1  Training  Accomplishments 

At  the  time  of  this  report,  as  the  recipient  of  the  Predoctoral  Traineeship  Award,  I  have  taken  20  out 
of  20  required  courses  towards  my  Ph.D.  degree  in  Medical  Physics.  The  courses  include  physics 
of  medical  imaging,  physics  of  radiation  therapy,  mathematics  for  medical  physicists,  statistics, 
anatomy  of  the  body,  radiation  biology,  and  teaching  assistant  training. 


2  Research  Accomplishments 

2.1  Investigation  of  non-circular  imaging  configurations 

Current  prototype  breast  CT  scanners  employ  a  circular  scanning  configuration  for  data  acquisi¬ 
tion  and  the  FDK  algorithm  for  approximate  image  reconstruction.  The  circular  scanning  configu¬ 
ration  is  easy  to  implement  and  can  provide  an  efficient  coverage  of  the  breast  when  a  flat-panel 
detector  is  used.  The  FDK  algorithm,  although  it  is  a  mathematically  approximate  algorithm,  re¬ 
mains  the  widely  used  algorithm  for  commercial  scanners  because  of  its  easy  implementation  and 
computational  efficiency.  However,  the  circular  imaging  configuration  yield  data  that  are  insuffi¬ 
cient  for  mathematically  exact  reconstruction  of  a  three-dimensional  (3D)  breast.  For  example, 
in  an  FDK  reconstruction,  images  within  two-dimensional  (2D)  slices  away  from  the  source  plane 
suffer  from  artifacts  such  as  contrast  reduction  and  streaks,  which  can  lead  to  an  increased  false 
positive  and/or  false  negative  diagnosis.  In  order  to  acquire  data  sufficient  for  accurate  3D  re¬ 
construction,  non-circular  scanning  configurations  have  been  proposed  and  investigated.  In  the 
first-year  research,  I  have  investigated  a  number  of  imaging  configurations,  as  described  below, 
that  can  yield  data  sufficient  for  mathematically  exact  image  reconstruction  of  the  breast. 

Helical  imaging  configuration:  The  helical  imaging  configuration  is  a  natural  extension  of 
the  circular  imaging  configuration  and  it  can  yield  data  sufficient  for  accurate  3D  reconstruction  of 
the  breast.  In  a  helical  configuration,  the  X-ray  source  and  detector  are  translated  while  they  are 
rotating,  and  the  source  trajectory  can  be  expressed  as 
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where  A  is  the  view  angle,  R  the  distance  from  the  source  point  to  the  rotation  axis,  and  h  the 
pitch  of  the  helical  trajectory,  which  is  defined  as  the  translation  distance  during  one  turn  of  rota¬ 
tion.  I  have  investigated  the  helical  configuration  using  both  simulated  and  real  data.  The  imaging 
configuration  of  our  simulation  study  is  identical  to  that  of  the  prototype  scanner  at  the  University 
of  California  ,  Davis.  In  the  configuration,  which  has  been  used  to  acquire  real  data  of  physical 
phantom  and  patient  breast,  the  distances  of  the  source  to  the  detector  plane  and  to  the  center 
of  rotation  are  87.78  cm  and  45.83  cm.  The  size  of  the  flat-panel  detector  is  39.7x29.7  cm2 
consisting  of  1024x768  detector  channels.  In  the  simulation  study,  I  have  used  different  digital 
phantoms,  such  as  the  Defrise  disk  phantom,  to  generate  data  at  500  views  using  analytic  pro¬ 
jection  equations  of  the  phantoms  for  every  turn  with  helical  pitches  ranging  from  3  cm  to  1 1  cm. 
I  used  the  disk  phantom  in  the  simulation  study  because  it  is  a  challenging  phantom  that  has 
been  used  widely  for  evaluation  of  image  reconstruction  algorithms.  I  investigated  the  volume 


coverage  that  can  be  reconstructed  exactly  using  different  helical  pitches.  The  simulation  study 
results  show  that  smaller  pitch  needs  more  turns  to  cover  the  same  volume  and  will  thus  result  in 
elevated  radiation  dose  to  the  patient  (assuming  dose  per  view  is  fixed).  However,  a  smaller  pitch 
can  result  in  better  coverage  for  the  volume  than  larger  helical  pitch. 

I  have  also  investigated  and  evaluated  the  helical  configurations  on  the  prototype  scanner. 
From  a  physical  Defrise  disk  phantom  and  a  sacrificed  mouse,  data  were  collected  by  use  of  four 
different  helical  pitches  ranging  from  6  cm  to  1 0  cm.  The  data  acquired  from  the  prototype  system 
requires  several  steps  of  processing  so  that  they  can  be  used  for  image  reconstruction.  Data 
processing  was  carried  out  for  removal  of  data  artifacts,  including  detector-response  correction, 
reference  projection  estimation,  and  physical  factor  correction.  The  detector-response  correction 
compensates  for  the  non-uniform  response  of  detector  channels  and  removes  dead/defect  de¬ 
tector  channels  if  correction  is  not  possible.  The  reference  projection  is  a  projection  acquired  in 
the  absence  of  the  imaged  object  and  is  required  for  the  conversion  of  the  raw  measurements  to 
data  for  image  reconstruction.  I  created  the  data  for  image  reconstruction  by  taking  the  logarithm 
of  the  ratio  between  reference  projection  and  object  projections.  I  have  estimated  the  reference 
projection  by  choosing  the  average  of  pixel  values  over  a  region  unattenuated  by  the  object  and 
using  the  value  as  a  uniform  reference  projection.  The  assumptions  behind  this  are  that  the  de¬ 
tector  response  is  uniform  and  that  the  object  is  sufficiently  small  so  that  its  projection  will  not  fully 
cover  the  detector.  I  have  also  compared  the  volume  coverage  in  real  data  studies  for  different 
helical  pitches.  In  order  to  cover  the  breast  volume  of  a  patient,  it  may  require  small  pitch  with 
multiple  turns  because  of  its  better  volume  coverage,  thus  radiation  dose  can  be  a  concern  for 
such  helical  configurations.  In  order  to  reduce  radiation  dose,  the  dose  at  each  view  for  a  helical 
configuration  need  to  be  reduced.  Given  the  fact  that  the  dose  level  per  view  is  already  low  and 
that  the  projection  data  have  significant  quantum  noise,  it  seems  unrealistic  to  further  reduce  the 
dose  for  each  view  in  a  helical  configuration. 

Saddle  imaging  configuration:  Saddle  imaging  configuration  is  another  configuration  that 
can  generate  data  sufficient  for  exact  3D  reconstruction  of  the  breast.  The  saddle  configuration 
can  be  interpreted  as  a  generalization  of  the  circular  configuration,  in  which  the  z-component  of 
the  X-ray  source  will  vary  periodically.  The  saddle  trajectory  under  study  can  be  expressed  as 

rb( A)  =  (i?cos  A,  i?sin  A, /icos2A)t  .  (2) 

I  have  studied  the  saddle  imaging  configurations  using  both  simulation  and  real  data.  The  geo¬ 
metric  parameters  of  the  simulation  system  is  the  same  as  what  was  used  in  the  helical  imaging 
configurations  except  for  that  different  values  of  h  were  used.  Again,  I  have  used  the  Defrise  disk 
phantom  to  generate  projection  data  at  500  view  over  2n  using  saddle  imaging  configurations.  I 
also  added  Gaussian  random  noise  terms  to  the  source  positions  in  order  to  simulate  mechanical 
inaccuracy  of  the  real  system.  I  have  compared  projection  data  and  reconstructed  images  for 
different  saddle  imaging  configurations. 

I  then  investigated  the  saddle  imaging  configuration  on  the  prototype  breast  CT  scanner.  The 
implementation  of  a  saddle  configuration  on  the  real  system  requires  the  X-ray  source  movement 
along  the  z-axis  according  to  a  non-linear  function  of  the  rotation  angle.  Because  of  mechanical 
limitations  of  motor  control,  such  a  non-linear  movement  can  not  be  implemented  directly  and  has 
to  be  realized  through  a  series  of  analytically  pre-determined  source  positions.  Additional  source 
positions  can  be  interpolated  between  two  adjacent  pre-determined  source  positions.  I  have 
computed  100  and  250  source  positions  distributed  equally  over  2n  using  an  analytic  formula  of 
the  trajectory.  Based  on  these  pre-calculated  source  positions,  I  interpolated  additional  400  or 
250  source  positions  for  additional  data  collection.  These  source  positions  were  used  for  the 
control  of  the  motor  and  data  have  been  collected  from  a  Defrise  disk  phantom  and  a  sacrificed 


mouse.  I  also  performed  detector  response  correction  and  reference  image  conversion  on  the 
projection  data,  and  have  investigated  the  volume  coverage  of  saddle  configurations.  The  results 
of  the  studies  indicate  that  a  saddle  configuration  may  have  a  good  coverage  for  the  plane  at 
2  =  0,  which  is  |  away  from  the  chest  wall. 

Circle-plus-line  imaging  configuration:  The  circle-plus-line  imaging  configuration  is  another 
configuration  that  can  yield  data  sufficient  for  exact  3D  reconstruction  of  the  breast.  It  can  be  re¬ 
alized  by  a  linear  scan  after  a  circular  scan,  and  it  can  be  readily  implemented  on  the  current 
prototype  scanners.  I  have  studied  this  imaging  configuration  using  simulation  data.  In  order  to 
cover  the  breast,  the  length  for  the  line  scan  can  be  up  to  20  cm.  The  configuration  requires  the 
full  coverage  of  the  breast  for  each  view,  including  the  linear  part.  Thus  requiring  large  detector 
size.  I  generated  data  at  500  views  for  the  circular  scan  and  500  views  for  the  line  scan  from  the 
digital  Defrise  disk  phantom,  and  compared  its  volume  coverage  to  that  of  the  helical  imaging  con¬ 
figuration  and  saddle  configuration.  Our  collaborators  have  also  implemented  the  circle-plus-line 
imaging  configuration  on  the  prototype  scanner,  in  which  the  Defrise  disk  phantom  was  scanned 
at  500  views.  I  have  performed  detector  response  correction  and  reference  image  conversion  on 
the  projection  data. 

Development  of  targeted  region-of-interest  imaging  approach:  I  have  conducted  a  pre¬ 
liminary  investigation  on  dynamic  collimation  to  achieve  targeted  ROI  imaging  using  simulation 
studies  [5],  in  which  the  source  and  the  detector  are  moved  for  each  view,  so  that  after  collima¬ 
tion,  the  ROI  can  always  be  covered  by  the  beam.  This  imaging  approach  leads  to  an  imaging 
problem  with  a  relatively  complex  source  trajectory  and  data  truncation,  which  can  however  be 
solved  with  the  back-projection  filtration  (BPF)  algorithm  [6,7],  I  have  used  the  BPF  algorithm  to 
reconstruct  ROI  images  from  truncated  data  collected  with  this  imaging  configuration. 

2.2  Investigation  and  development  of  novel  algorithms  tailored  to  accurate 
reconstruction  of  breast  CT  images 

Data  sufficiency  condition  for  breast  CT:  In  breast  CT,  a  half-cone  geometry  is  used  for  reduc¬ 
ing  the  radiation  dose  to  the  patient  torso.  This  can  significantly  alter  the  data  sufficiency  condition 
for  image  reconstruction  of  the  breast.  I  have  studied  the  data  sufficiency  conditions  using  helical 
and  other  imaging  configurations  for  both  full-cone  and  half-cone  geometries.  For  the  study  of 
half-cone  geometry,  I  first  generated  simulation  data  with  different  helical  pitches  from  the  Defrise 
disk  phantom  and  compared  the  reconstruction  volume  for  different  imaging  configurations.  When 
the  source  is  close  to  the  chest  wall,  I  removed  the  part  of  the  projection  data  that  will  be  blocked 
by  the  patient  body.  This  reduces  the  volume  that  can  be  reconstructed,  especially  the  part  near 
the  chest  wall.  The  study  results  suggest  that  the  use  of  a  helical  imaging  configuration  with  a 
half-cone  geometry  make  it  difficult  to  acquire  data  for  reconstructing  exactly  the  region  near  the 
chest  wall  of  the  patient. 

Noise  studies  of  image  reconstruction:  In  breast  CT  imaging,  because  of  the  radiation  con¬ 
cern,  the  total  permissible  radiation  dose  should  not  exceed  that  of  a  conventional  mammographic 
examination.  This  can,  however,  result  in  data  with  a  low  signal-to-noise  ratio  (SNR)  when  projec¬ 
tion  data  are  acquired  at  a  large  number  of  views.  In  addition  to  statistical  noise  due  to  low  dose, 
noise  is  also  introduced  in  the  signal  amplification  process  and  detector  electronics,  which  can 
significantly  impact  the  breast  CT  image  quality.  In  order  to  study  the  effect  of  noise,  I  have  gen¬ 
erated  noisy  projection  data  for  different  imaging  configurations,  with  different  levels  of  Gaussian 
noise  varying  from  0.5%  to  5%  to  represent  high  to  low  radiation  dose  levels.  I  have  investigated 
noise  properties  in  images  reconstructed  by  use  of  FDK  and  BPF  algorithms  at  different  noise 
levels  . 


2.3  Investigation  and  development  of  total-variation-based  algorithms  for 
breast  CT  image  reconstruction 

Radiation  dose  in  breast  CT  imaging  is  an  important  concern.  Further  dose  reduction  in  breast 
CT  can  be  achieved  through  collecting  data  at  a  reduced  number  of  projection  views.  When 
analytic  algorithms  such  as  the  FDK  and  BPF  algorithms  are  applied  to  sparse-view  data,  the 
reconstruction  images  will  contain  artifacts  such  as  streak  artifacts  that  can  lower  significantly 
the  image  utility.  Recently,  algorithms  based  upon  the  minimization  of  image  total  variation  (TV) 
subject  to  data  condition  have  been  developed  for  image  reconstruction  from  sparse-view  data 
[8,9].  I  have  conducted  a  preliminary  study  in  which  TV  algorithms  have  been  exploited  for  image 
reconstruction  from  sparse-view  data  in  breast  CT. 

Data  calibration  and  processing:  In  the  case  of  sparse-view  data,  data  information  avail¬ 
able  to  image  reconstruction  is  limited,  and  thus  it  is  critically  important  to  make  sure  that  data 
information  is  used  correctly.  The  first  issue  is  to  make  sure  that  data  geometry  is  calibrated  ap¬ 
propriately.  A  calibration  method  has  been  developed  to  obtain  data  geometric  parameters  from 
multiple  projections  of  metal  ball-bearing  (BB)  objects.  Assuming  that  the  system  has  a  mechan¬ 
ically  stable  rotation  center  and  that  the  detector  is  generally  within  the  rotation  plane,  one  can 
derive  the  geometrical  relationship  between  the  orbital  paths  of  the  individual  balls  for  determining 
data  geometric  parameters.  The  calibration  process  consists  of  identifying  the  center  position  of 
BBs  projection  at  each  view,  determining  the  ball  orbits,  and  relating  them  to  the  system  geome¬ 
try  parameters.  I  have  further  refined  and  validated  the  calibration  parameters  using  programs  for 
data  consistency  check  and  image  reconstruction.  Under  certain  conditions,  the  programs  can 
determine  the  projection  of  the  source  on  the  detector  plane  to  be  within  the  size  of  a  detector 
channel. 

TV-based  image  reconstruction:  Using  the  real  patient  breast  data  obtained  from  our  col¬ 
laborators  at  the  University  of  California  at  Davis,  I  have  performed  a  preliminary  investigation 
of  image  reconstruction  from  sparse-view  data  in  breast  CT  [10]  (Appendix  B).  In  the  study,  I 
have  added  projections  generated  from  small  spheres  that  simulate  micro-calcification  clusters 
(MCCs)  to  the  patient  data.  The  contrast  levels  of  the  simulated  MCC  structures  are  about  0.5 
to  1  times  of  that  of  the  background  glandular  tissues.  Furthermore,  I  have  investigated  noise 
properties  in  images  by  reconstructing  them  from  data  on  individual  detector  rows  or  from  data 
averaged  over  a  number  of  detector  rows.  In  the  latter  case,  the  noise  level  in  projection  data 
is  effectively  reduced.  The  patient  data  sets  contain  projections  collected  at  500  views  uniformly 
distributed  over  2ir.  In  the  sparse-view  data  study,  I  first  formed  sparse-data  sets  by  extracting 
from  500-view  patient  data  set  at  25  and  50  views  uniformly  distributed  over  2n.  I  then  used  both 
TV  and  FDK  algorithms  to  reconstruct  images  from  the  extracted  sparse-view  data  sets.  The 
reconstruction  results  indicate  that  the  TV  algorithm  can  yield  generally  higher  contrast  levels  of 
the  MCC  structures  than  does  the  FDK  algorithm.  This  conclusion  applies  to  both  data  with  and 
without  averaging  as  discussed  above.  An  implication  of  the  study  is  that  the  TV  algorithm  can  be 
employed  for  potentially  reducing  substantially  the  total  radiation  dose  in  breast  CT. 

2.4  Evaluation  of  the  proposed  configurations  and  algorithms 

The  evaluation  study  is  planned  largely  for  years  2  and  3.  However,  I  have  carried  out  preliminary 
studies  on  data  and  image  noise  properties  this  year  to  lay  out  the  basis  for  the  planned  evaluation 
studies  in  the  next  year. 

Digital  breast  phantom:  One  of  the  important  tools  in  evaluation  studies  is  the  availability  of 
digital  breast  phantoms  containing  structures  with  different  levels  of  complexity.  I  have  created 


and  modified  two  breast  phantoms  that  consist  of  a  half-ellipsoid  object  wrapped  by  simulated 
skin.  The  half-ellipsoid  contains  a  compound  of  adipose  and  glandular  tissue,  which  is  also  ap¬ 
proximated  by  ellipsoids  of  different  sizes.  Spheres  with  strong  attenuation  (i.e.,  high  contrast) 
of  different  sizes  that  mimic  MCCs  have  also  been  included  in  the  phantoms.  Specifically,  in  a 
typical  MCC,  the  smallest  calcification  is  about  200  n m,  and  a  MCC  typically  includes  3  to  10 
simulated  micro-calcifications.  The  MCCs  have  been  added  to  the  glandular  tissues  and  the  axil¬ 
lary  tail  of  the  breast.  From  the  digital  breast  phantoms,  I  have  generated  projection  data  at  500 
views  uniformly  distributed  over  27t  using  a  circular  imaging  configuration  that  is  currently  used  by 
our  collaborators.  I  have  also  generated  noisy  data  by  adding  Gaussian  noise  to  the  simulated 
data.  The  digital  breast  phantoms  and  simulation  tools  developed  provide  a  solid  basis  for  the 
evaluation  study  planned  for  the  next  year. 

Noise  study  of  ROI-imaging  approach:  I  have  studied  the  noise  properties  in  images  recon¬ 
structed  by  use  of  the  BPF  algorithm  in  a  simulation  study,  in  which  a  uniform  cylinder  phantom  is 
used  to  generate  1 000  noisy  data  sets  containing  Gaussian  noise.  From  the  1 000  noisy  data  sets, 
I  first  reconstructed  noisy  images  using  the  BPF  algorithm  and  then  computed  empirical  image 
variances  and  covariances.  The  results  show  that  the  noise  level  within  an  ROI  is  non-uniform 
and  that  the  image  noise  level  in  the  center-of-rotation  region  is  generally  lower  than  that  in  other 
regions.  Therefore,  it  is  important  to  design  an  imaging  configuration  such  that  the  ROI  can  be 
placed  near  the  center  of  rotation. 


KEY  RESEARCH  ACCOMPLISHMENTS 


•  I  have  investigated  and  evaluated  helical,  saddle,  and  circle-plus-line  imaging  configurations 
for  potential  use  in  breast  CT. 

•  I  have  carried  out  a  preliminary  investigation  on  dynamic  collimation  for  ROI  imaging  using 
simulation  data. 

•  I  have  modified  and  developed  a  BPF  algorithm  for  image  reconstruction  from  data  acquired 
with  non-circular  imaging  configurations  in  breast  CT. 

•  I  have  modified  and  developed  a  BPF  algorithm  for  image  reconstruction  from  simulation 
data  acquired  in  an  imaging  configuration  employing  dynamic  collimation. 

•  In  working  with  the  collaborators,  projection  data  of  physical  phantoms  have  been  collected. 

•  In  working  with  the  collaborators,  projection  data  of  patient  breasts  have  been  collected. 

•  I  have  performed  a  preliminary  investigation  on  potential  dose  reduction  by  lowering  the 
number  of  projection  views  in  breast  CT. 

•  I  have  created  digital  breast  phantoms  with  different  structures  and  micro-calcification  clus¬ 
ters,  which  can  be  used  for  quantitative  evaluation  of  imaging  configurations  and  reconstruc¬ 
tion  algorithms. 
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CONCLUSIONS 


As  the  recipient  of  the  Predoctoral  Traineeship  Award,  I  have  finished  the  required  courses 
towards  my  Ph.D.  degree.  These  trainings  have  proven  useful  for  me  to  achieve  the  proposed 
research  goals. 

During  the  first  year,  I  have  investigated  and  evaluated  non-circular  imaging  configurations 
potentially  useful  for  breast  CT  imaging  by  studying  data  sufficiency  conditions  and  by  performing 
both  simulation  and  real  data  studies.  I  have  developed  and  modified  BPF  algorithms  to  recon¬ 
struct  images  within  an  ROI  from  data  containing  truncations.  The  modified  BPF  algorithm  can 
yield  more  accurate  images  than  the  existing  FDK  algorithm.  In  an  attempt  to  further  reduce 
imaging  dose,  one  may  reduce  the  number  of  views  at  which  projection  data  are  collected.  In  this 
situation,  analytic  algorithms  may  not  provide  adequate  results.  Instead,  I  have  developed,  modi¬ 
fied,  and  implemented  the  TV  algorithm  to  reconstruct  images  from  sparse-viewdata.  Also,  I  have 
created  digital  breast  phantoms  that  can  be  used  as  an  important  tool  for  evaluation  of  the  poten¬ 
tial  utility  of  imaging  configurations  and  reconstruction  algorithms.  Furthermore,  I  have  performed 
a  preliminary  evaluation  study  of  different  reconstruction  algorithms  in  which  both  simulation  and 
real  data  were  used. 

In  summary,  I  believe  that  I  have  achieved  the  goals  planned  for  the  first  year  and  laid  down  the 
foundation  for  the  research  in  the  next  two  years.  The  aims  in  the  next  two  years  include  further 
development  of  the  TV  algorithm  for  image  reconstruction  in  low-dose  breast  CT,  investigation 
of  the  effects  of  physical  factors  on  the  quality  of  image  reconstruction,  and  performance  of  well 
designed  evaluation  studies. 
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Abstract:  The  back-projection-filtration  (BPF)  algorithm  has  been  applied  to  image  reconstruction  for 
cone-beam  configurations  with  general  source  trajectories.  The  BPF  algorithm  can  reconstruct  3-D  region- 
of-interest  (ROI)  images  from  data  containing  truncations.  However,  like  many  other  existing  algorithms  for 
cone-beam  configurations,  the  BPF  algorithm  involves  a  back-projection  with  a  spatially  varying  weighting 
factor,  which  can  result  in  the  non-uniform  noise  levels  in  reconstructed  images  and  increased  computation 
time.  In  this  work,  we  propose  a  BPF  algorithm  to  eliminate  the  spatially  varying  weighting  factor  by  using  a 
rebinned  geometry  for  a  general  scanning  trajectory.  This  proposed  BPF  algorithm  has  an  improved  noise 
property,  while  retaining  the  advantages  of  the  original  BPF  algorithm  such  as  minimum  data  requirement. 

Keywords:  computed  tomography;  region  of  interest;  general  trajectory;  rebinned  backprojection-filtration 


Introduction 

Significant  advances  have  been  made  in  the  develop¬ 
ment  of  theoretically  exact  algorithms  for  image  re¬ 
construction  from  cone-beam  projections.  Theoretically 
exact  algorithms  that  were  developed  initially  for  a 
helical  trajectory [1  "3]  have  been  extended  to  reconstruct 
images  from  cone-beam  projections  acquired  with 
general  source  trajectories14"81.  The  back-projection- 
filtration  (BPF)  algorithm  is  one  of  the  recently  devel¬ 
oped  algorithms19"1 1].  It  can  be  applied  to  reconstruct¬ 
ing  images  for  a  wide  class  of  general  scanning  trajec¬ 
tories  and  is  capable  of  reconstructing  an  image  within 
a  region-of-interest  (ROI)  from  projection  data  con¬ 
taining  truncations.  This  property  of  the  BPF  algorithm 
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allows  imaging  applications  of  practical  significance. 
For  example,  in  non-conventional  computed  tomogra¬ 
phy  (CT)  applications,  it  is  not  uncommon  that  the 
available  detector  covers  only  a  portion  of  a  field  of 
view  (FOV)  that  is  needed  otherwise  for  completely 
covering  the  entire  support  of  the  imaged  object.  More¬ 
over,  data  acquisition  in  many  practical  applications  is 
achieved  through  the  rotation  of  the  object  around  the 
physical  center  of  rotation,  which  is  often  chosen  as  the 
center  of  mass  of  the  object  and  may  be  at  a  distance 
from  the  FOV  center,  as  shown  in  Fig.  1.  Although  this 
imaging  approach  leads  to  an  imaging  problem  with  a 
relatively  complex  source  trajectory  and  data  truncation, 
one  can  apply  the  BPF  algorithm  to  reconstructing  an 
ROI  image  from  truncated  data  collected  with  this  kind 
of  imaging  configurations  with  general  trajectories1121. 

Like  many  existing  algorithms  for  divergent-beam 
configurations,  the  original  BPF  algorithm  also  in¬ 
volves  the  computation  of  a  spatially  varying  weight¬ 
ing  factor  in  its  back-projection  step,  which  can  result 
in  non-uniform  noise  levels  in  reconstructed  images 
and  increased  computation  time113"161.  Therefore,  it  is 
desirable  to  eliminate  the  spatially  varying  weighting 
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Fig.  1  Illustration  of  the  3-D  imaging  configuration. 
The  object  rotates  around  the  rotation  axis  O.  The 
source  can  move  freely  vertically  and  horizontally 
while  the  detector  can  only  move  horizontally.  The  re¬ 
gion  enclosed  by  the  dashed  curve  indicates  the  FOV 
covered  by  the  detector. 

factor  in  the  BPF  algorithm  for  improving  its  noise 
properties  and  computational  efficiency.  In  this  work, 
based  upon  the  original  BPF  algorithm,  we  propose  a 
BPF  algorithm  in  terms  of  a  rebinned  geometry  for  a 
general  scanning  trajectory  designed  for  ROI  imaging, 
in  which  no  spatially  varying  weighting  factor  is  in¬ 
volved  in  the  back-projection  step.  The  use  of  the  re¬ 
binned  geometry  can  not  only  eliminate  the  spatially 
varying  weighting  factor,  thus  improving  the  noise 
properties  of  BPF  reconstructions,  but  also  retain  the 
properties  of  the  original  BPF  algorithm  such  as 
minimum  data  requirement  and  ROI-image  reconstruc¬ 
tion  from  truncated  data. 

1  New  ROI  Imaging  Approach 

1.1  Design  of  an  ROI  imaging  approach 

We  display  in  Fig.l  a  3 -dimensional  (3-D)  scanning 
configuration  under  consideration,  in  which  the  object 
rotates  around  the  rotation  axis,  the  source  can 
move  freely  vertically  and  horizontally,  and  the  detec¬ 
tor  can  only  move  horizontally.  The  line  connecting  the 
source  and  the  detector  center  is  always  perpendicular 
to  the  detector  plane.  A  sketch  of  this  scanning  geome¬ 
try  within  2-dimensional  (2-D)  slice  at  z  =  0  is 
shown  in  Fig.  2.  An  ellipsoid  support  enclosed  by  the 
thin  curve  represents  the  imaged  object  in  which  an 
ROI  is  indicated  as  the  shaded  region,  and  the  FOV  is 


enclosed  by  a  thick  curve.  The  geometric  center  of  the 
ellipsoid  is  placed  at  the  physical  center  of  rotation  of 
the  imaging  system,  which  is  also  chosen  as  the  origin 
of  the  fixed  coordinate  system.  We  assume  that  the 
distance  of  the  FOV  center  to  the  physical  center  of 
rotation  is  R0,  that  the  width  of  the  detector  is  2D, 

and  that  the  distance  H  between  the  physical  center 
of  rotation  and  the  detector  remains  constant.  Let  R 
denote  the  radius  of  the  FOV.  Clearly,  R<D  for  the 
cone -beam  projection  under  consideration.  Projections 
at  different  views  are  collected  through  the  rotation  of 
the  object  about  the  physical  rotation  center  and  the 
motion  of  the  source  and  detector  in  the  space,  which 
is  constrained  by  two  conditions:  (1)  the  size  of  the 
FOV  (enclosed  by  the  thick  curve)  formed  by  the 
source  and  detector  remains  unchanged  and  always 
covers  the  ROI,  and  (2)  the  source  point,  the  FOV 
center,  and  the  detector  mid-point  always  remain  on 
the  same  line,  which,  without  loss  of  generality,  is  as¬ 
sumed  to  be  perpendicular  to  the  detector  plane.  An 
additional  constraint  on  the  source  is  that  it  is  outside 
the  rotating  object. 


Source 


Fig.  2  Illustration  of  the  scanning  geometry  within 
2-D  slice  at  z  =  0  .  The  support  of  the  imaged  object  is 
enclosed  by  the  thin  curve.  During  imaging,  the  object 
rotates  around  the  physical  center  of  rotation.  The  ROI 
(shaded  area)  to  be  imaged  is  within  the  FOV  (enclosed 
by  the  thick  curve)  formed  by  the  source  and  the  de¬ 
tector.  When  the  object  is  rotated,  the  detector  can 
move  only  along  x-axis,  whereas  the  source  can  move 
along  x-  and/or  y-axes  so  that  the  FOV  always  cov¬ 
ers  the  ROI  completely. 

1.2  Source  trajectory 

When  the  imaged  object  rotates,  the  source  and  the 
detector  are  moved  accordingly  so  that  they  form  an 
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FOV  that  always  covers  the  ROI.  Consequently,  the 
location  of  the  FOV  center  also  changes.  Considering 
the  motion  constraints  on  the  source  and  detector  de¬ 
scribed  above,  one  can  determine  the  trajectories  of  the 
source,  the  FOV  center,  and  the  mid-point  of  the  de¬ 
tector  in  the  x-y  plane.  Let  A  denote  the  rotation  an¬ 
gle  of  the  object,  which  is  defined  as  the  angle  between 
the  long  axis  of  the  ellipsoid  object  and  the  x  axis. 
Therefore,  the  trajectories  of  the  FOV  center 
(xF,yF,zF)  and  detector  center  ( xD,yD,zD )  in  the 


fixed  coordinate  system  can  be  expressed  as 

(xF ,  y¥ ,  zF )  =  (R0  co  s  X,  R0  sin  2, 0) 

(xD,yD,zD)  —  (7?q  cos  2,  —  H ,  0) 

Also,  one  can  express  the  source  trajectory 

the  fixed  coordinate  system  as 

r,W  =  (xs,ys,zs)  =  (R0  cosA,ys0,0) 

where  vs0  =  D tan(#S,  +  </>,)- H,  and 

,  H  +  Rn  sin  A 

m,  =  arctan - - - 

1  D 

,  R 

<p2  =  arctan— p 

^D2 +{H +  R0smA)2 -R2 


(1) 

(2) 

rs(X)  in 

(3) 

(4) 

(5) 


As  shown  in  Fig.  2,  ^  is  the  angle  between  the  x  axis 


and  the  line  connecting  the  FOV  center  and  the  detec¬ 
tor  edge,  and  </>2  indicates  the  angle  between  the  line 


connecting  the  FOV  center  and  the  detector  edge  and 
the  line  connecting  the  source  and  the  detector  edge. 

Considering  the  physical  constraint  that  the  source 
cannot  be  within  the  object,  we  assume  that  the  object 
support  has  a  cylindrical  shape,  which  has  the  same 
middle  transverse  slice  as  the  ellipsoid  object.  Assum¬ 
ing  that  the  long  axis  of  the  ellipsoid  object  is  aligned 
with  the  x  axis  when  X  =  0,  we  can  rewrite  the  source 

trajectory  in  Eq.  (3)  as 

rs(X)  =  (xs,ys,zs)  = 

(R0cosX,R0sinX  +  max(ys0  - R0sinX,t),0)  (6) 
where 

t_-b2R0smX  +  T(a,b,R0-,X) 
a2  cos2  X+b2  sin2  X 


and 

T (a, b,  R0;X)=^j  a4  b2  cos2  X-a2b2Rg  cos2  X+a2b4  sin2  X 

(8) 

In  general,  reconstruction  algorithms  are  developed 
for  imaging  configurations  in  which  the  object  is  fixed. 
Therefore,  we  need  to  determine  the  source  trajectory 
in  a  coordinate  system  in  which  the  object  is  fixed  so 


that  these  algorithms  can  be  applied  directly.  In  the 
case  under  study,  we  refer  to  the  coordinate  system 
fixed  on  the  object  as  a  object- fixed  coordinate  system 
(jc0,  v„,z0) ,  which  can  be  related  to  the  original  coor¬ 
dinate  system  (x,y,z)  described  above  as 

cosA  -sinA  0 

Oo>  To>  zo)  =  (v  y>  z ) 


sin  A 
0 


cos  A 
0 


(9) 


Using  jcs,  ys,and  zs  to  replace  x,  y,  and  z  in  the 
right-hand-side  of  Eq.  (9),  one  can  obtain  the  expres¬ 
sion  of  the  source  trajectory  ro(X)  =  (x0,y0,z0)  in  the 
object- fixed  coordinate  system,  which  is  displayed  in 
Fig.  3. 


Fig.  3  Illustration  of  the  scanning  geometry  in  the 
object-fixed  coordinate  system.  The  thin  solid  and 
dashed  lines  represent  the  source  trajectory  and  de¬ 
tector  center  trajectory  respectively.  The  FOV  is  en¬ 
closed  by  the  thick  curve  (shaded  area). 


2  BPF  Algorithm  for  ROI-Image 
Reconstruction 


2.1  Original  BPF  algorithm 

In  order  to  describe  the  original  BPF  algorithm,  we 
introduce  a  rotation-coordinate  system  {m,v,w}> 

whose  origin  is  fixed  on  the  source  point  for  the  con¬ 
venience  of  the  reconstruction.  Its  three  unit  vectors  in 
the  fixed-coordinate  system  are 

eu  =  (cosA,  -sin A,  0)T, 

e„  =  (0,  0,  1)T, 

e„  =  (sin  A,  cos  A,  0)T  (10) 

Consider  a  flat-panel  detector  with  its  normal  direc¬ 
tion  along  ew  and  at  a  distance  S  from  the  source 
point.  A  2-D  coordinate  system  {u,v}  is  assumed  to 
be  fixed  on  the  detector  plane  and  the  u  and  v  axes  are 
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along  the  unit  vectors  eu  and  £v .  Any  point  on  the 
detector  can  thus  be  specified  by  two  parameters  u 
and  v.  At  source  position  A,  the  cone-beam  projec¬ 
tion  of  the  object  function  f(r)  at  a  point  («d,vd) 
on  the  detector  can  be  expressed  as 

P(ui,vi,A)  =  \of(ro(A)  +  tj3)dt  (11) 


where 

P  =  I  ,  \  ,Kg.W  +  vdev(A)  -  SeJA)\  (12) 

M+vl+s 

is  a  unit  vector  indicating  the  direction  of  the  ray  that 
starts  from  source  points  roW  and  passes  through 
the  point  (w,v)  on  the  detector. 

The  original  BPF  algorithm  can  reconstruct  ROI 
images  from  truncated  data  acquired  with  a  general 
trajectory  by  reconstructing  images  on  chords.  For  a 
given  continuous  trajectory,  a  chord  is  a  line  segment 
comiecting  two  points  ;-0(Aa)  and  r 0(Tb)  on  the 


trajectory.  The  fonn  of  the  original  BPF  algorithm  used 
here  is  given  in  Ref.  [8],  and  the  image  on  the  chord 
specified  by  Ta  and  Ab  is  given  by 


f'<*"A"A>>  2tC  b(xc) ' 

^bOOKgM,Aa,Ab)  +  2nP0 


Xc  ~  Xc 

where  b(xc )  =  A J(xc2  -  xc)(xc  -  xcl)  and  P0  denotes 
the  projection  data  along  the  chord.  gc(xc,Aa,/ tb)  is 


defined  by  the  following  equation, 

f  3  3  x  rA  d A  d 

gAx  ,,/l.,/lh)  = - 

Ja*  Hxc)-r,{A)\AA 


P(ud,v6,A) 


(13 


) 

By  use  of  the  original  BPF  algorithm,  one  can  re¬ 
construct  ROI  images  from  the  truncated  data  acquired 
with  the  source  trajectory  described  in  Section  1.  It  can 

be  observed  that  the  weighting  factor  — — — — —— 

Ir0cc)_fo(/Ol 


in  the  back-projection  step  (i.e.,  Eq.  (13))  is  spatially 
varying,  which  can  result  in  the  increased  noise  level 
and  computation  load.  We  seek  to  eliminate  such  a 
factor  by  using  the  rebinned  geometry,  as  discussed 
below. 


2.2  Rebinned  BPF  algorithm 

In  the  rebinned  BPF  algorithm,  the  acquired  projection 


data  are  first  rebinned  into  the  fan-parallel-beam 
geometry  [17‘20].  The  rebinned  data  can  be  expressed  as 
Q'(uA,vA,(p)  satisfying 

Q'(ua,va,<p)  =  P(ua,va,A)  (14) 

given  that 

u 

cp  =  A-  arctan—  (15) 

S 

By  using  this  relationship,  the  weighting  factor 

— 7—7  in  Eq.  (13)  can  be  eliminated. 
|r(xc)-roWI 

3  Numerical  Studies 


We  have  performed  computer-simulation  studies  to 
investigate  and  evaluate  the  rebinned  BPF  algorithm 
for  achieving  an  ROI  reconstruction  for  the  trajectory 
described  in  Section  2.  The  numerical  phantom  is 
modified  from  a  standard  Shepp-Logan  phantom.  The 
center  of  the  standard  Shepp-Logan  phantom  is  first 
shifted  to  the  center  of  the  FOV  and  the  size  of  the 
Shepp-Logan  phantom  is  then  scaled  to  fit  the  size  of 
the  FOV.  A  long  narrow  ellipsoid  object  is  placed  at  the 
rotation  center  which  connects  the  Shepp-Logan  phan¬ 
tom  to  the  rotation  axis.  To  cover  the  ROI,  the  FOV 
has  a  radius  R  =30  cm  and  is  at  a  distance  R0  =  50  cm 

from  the  physical  center  of  rotation,  and  the  detector 
has  a  size  of  D  =  64  cm  and  is  placed  at  a  distance 
from  the  x  axis  at  77  =  100  cm.  The  detector  has 
512x512  units.  Using  these  parameters  in  the  imag¬ 
ing  configuration,  we  computed  cone-beam  projection 
data  from  the  rotating  object  at  720  views  from  -A 


to 
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from  which,  noisy  data  were  generated  by 


adding  Gaussian  noise  with  a  standard  deviation  that  is 
about  1.5%  of  the  maximum  projection  value.  From 
these  noiseless  and  noisy  data  sets,  we  have  recon¬ 
structed  the  ROI  images.  In  Fig.  4,  we  displayed  the 
images  reconstructed  from  the  noisy  data  by  use  of 
rebinned  and  original  BPF  algorithms. 

We  have  also  perfonned  a  study  to  investigate  the 
noise  properties  of  the  rebinned  BPF.  In  this  study,  a 
uniform  cylinder  phantom  was  used.  From  the  uniform 
phantom,  we  generated  1000  sets  of  the  data  contain¬ 
ing  stationary  Gaussian  noise.  Without  loss  of  general¬ 
ity,  we  focus  on  reconstructing  images  on  a  set  of 


chords  specified  by  Al  =  ^  and  A2  = 


5  71 
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the  1000  images  on  these  chords  reconstructed  from 
the  noisy  data  sets,  we  calculated  the  empirical  image 
variances  on  the  chords.  For  comparison,  the  empirical 
image  variances  are  obtained  from  the  reconstructed 
images  by  use  of  the  original  BPF  algorithm.  The  im¬ 
age  variances  along  the  chord  specified  by  \  ^ 


and  A2  =  are  displayed  in  Fig.  5.  It  can  be  ob¬ 


served  that  the  rebinned  BPF  algorithm  generally 
yields  images  with  lower  and  more  spatially  uniform 
variances  than  does  the  original  BPF  algorithm.  This 
uniform  noise  property  in  images  obtained  by  use  of 
the  rebinned  BPF  algorithm  is  a  direct  result  of 
the  elimination  of  the  spatially  varying  weighting 
factor  from  the  back-projection  in  the  rebinned  BPF 
algorithm. 


•  00 

©  ©  o 


Fig.  4  The  noisy  images  obtained  with  the  rebinned 
BPF  algorithm  (top  row)  and  the  original  BPF  algo¬ 
rithm  (bottom  row)  within  the  2-D  slices  at  a=0  mm, 
>’=()  mm,  and  -=0  mm.  The  display  window  is  [0.98, 
1.2], 


Fig.  5  Image  variances  on  a  chord,  specified  by 
ij=-^  and  i2  =  ^  ,  obtained  by  use  of  the  re¬ 
binned  BPF  algorithm  (dashed  line)  and  original  BPF 
algorithm  (solid  line) 


4  Conclusions 

In  the  work,  we  have  described  a  rebinned  BPF  algo¬ 
rithm,  which  involves  no  spatially  varying  weighting 
factor  in  its  back-projection  step,  for  ROI-image 
reconstruction  from  truncated  data  acquired  with  a 
general  trajectory.  We  have  performed  computation- 
simulation  studies  to  validate  and  evaluate  the  rebinned 
BPF  algorithm.  The  quantitative  results  demonstrate 
that  the  exact  ROI-image  reconstruction  can  be  ob¬ 
tained  with  the  rebinned  BPF  algorithm.  Most  impor¬ 
tant,  the  rebinned  BPF  algorithm  can  improve  the  noise 
properties  in  terms  of  image  variances.  It  may  find 
practical  implications  in  numerous  non-conventional 
CT  scans  involving  general  trajectories. 
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I  Introduction 

Breast  cancer  is  the  most  frequently  diagnosed  can¬ 
cer  and  the  second  leading  causes  of  cancer  mortality 
among  woman.  Researchers  are  actively  developing 
x-ray  breast  imaging  techniques  that  can  form  three- 
dimensional  image  of  the  breast,  including  breast  to¬ 
mosynthesis  and  dedicated  breast  computer  tomogra¬ 
phy  CT.  Several  prototype  dedicated  breast  CT  scan¬ 
ners  have  been  built  in  different  institutions  through¬ 
out  US.  Most  of  the  scannesr  have  a  circular  configura¬ 
tion  and  adopted  an  approximate  analytic  reconstruc¬ 
tion  algorithm,  which  requires  a  large  number  of  views 
to  be  acquired.  The  total  dose  of  a  breast  CT  scan  is 
normally  about  the  same  as  a  two- view  mammogram. 
It  is  still  desirable  to  further  reduce  the  breast  CT 
dose  for  general  screening  purpose. 

Inspired  by  the  principle  of  compressive  sensing, 
we  have  proposed  algorithms  for  image  reconstruc¬ 
tion  from  fan-  and  cone-beam  data  collected  at  highly 
sparse  views  through  minimization  of  the  total  varia¬ 
tion  (TV)  of  the  image  subject  to  the  condition  that 
the  estimated  data  is  consistent  to  the  measured  data. 
The  advantage  of  our  TV  minimization  algorithm  is 
that  by  reducing  the  number  of  views,  total  dose  can 
be  reduced  significantly.  In  this  study,  we  investigate 
the  sparse  data  reconstruction  for  breast  C'T  by  using 
very  few  projections.  The  aim  of  this  study  is  to  per¬ 
form  a  preliminary  investigation  to  further  reduce  the 
dose  of  breast  CT  for  breast  cancer  screening. 

II  Materials  and  Methods 

The  breast  C'T  data  were  acquired  by  a  clinical  trial 
conducted  by  our  collaborator  at  UC-Davis  medical 
center.  The  prototype  scanner  consists  of  a  water- 
cooled  tungsten  anode  x-ray  tube,  a  flat-panel  Csl  in¬ 
direct  x-ray  detector  that  was  40  cm  wide  and  30  cm 
tall  and  an  integrated  bearing-motor-encoder  system. 
The  detector  has  a  pixel  matrix  of  1024  x  768,  with 
pixel  size  of  0.388  mm.  The  source  to  detector  dis¬ 
tance  is  87.78  cm  and  the  source  to  rotation  center 
is  45.83  cm.  For  each  data  sets,  500  cone-beam  pro¬ 


jection  images  are  acquired  360°  around  the  patient’s 
breast  in  16.6  seconds. 

In  the  TV  algorithm,  the  imaging  system  is  mod¬ 
eled  by  a  discrete  linear  system  g  =  Af/,  where  g 
is  the  measured,  /  denotes  the  discrete  image  and 
M  is  the  system  matrix.  The  TV  algorithm  is  for¬ 
mulated  into  a  constrained  TV-minimization  frame¬ 
work,  which  seeks  to  find  image  /  that  minimizes  its 
TV:  f*  =  argmin  /  ,  subject  to  the  data  fidelity 

and  non-negativity  constraints.  The  proposed  TV  al¬ 
gorithm  solves  this  constrained  optimization  problem 
through  using  a  hybrid  scheme  in  which  the  projec¬ 
tion  onto  convex  sets  (POCS)  is  used  to  enforce  the 
data  constraint,  and  gradient  descent  method  is  used 
subsequently  for  minimizing  the  image  TV. 

In  order  to  investigate  the  performance  of  the  TV 
algorithm,  we  extract  the  mid-slice  sinogram  from  the 
patient  data  and  add  simulated  sinogram.  The  simu¬ 
lated  sinogram  is  generated  by  analytic  phantom  with 
exactly  the  same  geometry  as  the  prototype  breast  CT 
scanner.  The  analytic  phantom  consists  of  small  ellip¬ 
sis  to  mimic  small  structures  in  the  breast  with  their 
sizes  varying  from  0.5  mm  to  1.3  mm.  The  intensity  of 
the  structures  are  about  0.5  to  1.  times  of  the  intensity 
of  the  background  glandular  tissues.  The  modified  pa¬ 
tient  breast  image  is  shown  in  the  first  collumn  of  Figs 

1  to  4  and  the  white  dots  show  the  small  structures, 
the  intensity  of  the  small  structures  are  about  1.5  to 

2  times  of  the  background  tissues. 

In  order  to  study  the  effects  of  the  noise,  we  also 
average  mid-slice  sinogram  with  its  two  neighboring 
detector  rows  to  form  the  new  mid-slice  sinogram  and 
add  simulated  sinogram  to  this  new  sinogram.  By  av¬ 
eraging,  the  noise  level  in  the  projection  data  is  effec¬ 
tively  reduced.  From  the  full  500  views  of  projection 
data  sets  of  the  mid-slice  sinogram  with  and  without 
averaging,  we  extract  25  and  50  views  of  projection 
data  which  are  uniformly  distributed  from  0  to  360°. 
We  used  both  TV  and  FDK  algorithms  to  reconstruct 
images  from  these  four  sparse  data  sets.  The  recon¬ 
struction  image  size  is  15  x  15  cm2and  the  correspond¬ 
ing  image  array  size  is  512  x  512. 


1 


Because  the  total  dose  of  breast  CT  need  to  be  dis¬ 
tributed  over  a  large  number  of  views  (usually  several 
hundred),  there  are  considerable  quantum  noise  in  the 
images  reconstructed  directly  by  FDK  algorithms  .  So 
we  applied  a  hanning  filter  to  the  images  reconstructed 
by  FDK  algorithms.  For  TV,  we  take  the  images  at 
the  40th  iteration. 


Fig.  2.  Image  reconstructed  by  TV  (middle)  and  FDK  (right)  algorithms 
using  25  views  of  projection  data  from  the  averaged  sinogram.  The  mod¬ 
ified  breast  image  (left)  shows  the  small  structures.  The  display  window 
is  [0.15,0.28]  cm-1. 


Ill  Results 


In  Figs.  1  to  4,  we  displayed  the  images  recon¬ 
structed  by  use  of  TV  and  FDK  for  the  mid-slice  sino¬ 
gram  and  the  averaged  sinogram.  In  Fig  1.,  the  25- 
view  FDK  reconstruction  is  clearly  overwhelmed  by 
the  streaks  caused  by  aliasing  and  it  is  very  difficult 
to  differentiate  between  real  structures  and  artifacts. 
However,  the  image  reconstructed  by  TV  algorithm  is 
free  of  streak  artifacts  and  some  of  the  structures  can 
be  easily  identified  .  In  Fig  2.,  we  displayed  the  TV 
and  FDK  reconstruction  from  the  25-view  averaged 
sinogram.  Both  the  FDK  reconstruction  and  TV  re¬ 
construction  become  better  in  this  case.  However,  the 
FDK  image  is  still  dominated  by  the  streak  artifacts, 
while  almost  all  four  groups  of  the  small  structures  can 
be  identified  in  the  TV  reconstruction  image.  In  Fig 
3,  we  displayed  the  TV  and  FDK  reconstruction  from 
the  50-views  of  projection  data  without  averaging  . 
There  are  still  artifacts  caused  by  aliasing  and  the  ar¬ 
tifacts  have  very  close  appearance  to  the  small  struc¬ 
tures,  which  will  caused  many  false  positives.  While 
in  the  TV  reconstruction  images,  the  structures  can 
be  visualized  much  better  and  all  the  structures  can 
be  identified  in  the  image.  In  Fig.  4,  we  displayed  the 
TV  and  FDK  reconstruction  from  the  50-views  aver¬ 
aged  sinogram.  The  visualization  of  the  structures  are 
almost  the  same  in  this  case. 


Fig.  1.  Image  reconstructed  by  TV  (middle)  and  FDK  (right)  algo¬ 
rithms  using  25  views  of  projection  data  without  averaging.  The  modified 
breast  image  (left)  shows  the  small  structures.  The  display  window  is 
[0.15,0.28]  cm-1. 


Fig.  3.  Image  reconstructed  by  TV  (middle)  and  FDK  (right)  algo¬ 
rithms  using  50  views  of  projection  data  without  averaging.  The  modi¬ 
fied  breast  image  (left)  shows  thesmall  structures.  The  display  window  is 
[0.15,0.28]  cm-1. 


Fig.  4.  Image  reconstructed  by  TV  (middle)  and  FDK  (right)  algorithms 
using  50  views  of  projection  data  from  the  averaged  sinogram.  The  mod¬ 
ified  breast  image  (left)  shows  the  additional  small  structures.  The  dis¬ 
play  window  is  [0.15,0.28]  cm-1. 

IV  Discussion 

A  preliminary  study  of  sparse  data  reconstruction 
for  breast  CT  using  TV  minimization  algorithm  is  per¬ 
formed  by  using  real  patient  data  with  adding  small 
structures  added  to  the  breast.  By  comparing  the 
FDK  reconstruction  and  TV'  reconstruction,  it  can  be 
seen  that  the  structures  can  be  visualized  much  better 
in  image  reconstructed  by  TV  algorithm  than  the  im¬ 
ages  reconstructed  by  FDK  for  25  views  of  projection 
data  of  the  two  different  noise  levels  since  the  streak 
artifacts  dominated  the  FDK  reconstruction.  For  im¬ 
ages  reconstructed  from  50views  of  projection  data, 
the  streak  artifacts  are  not  a  serious  issue  any  more. 
However,  the  TV  reconstruction  is  still  better  in  vi¬ 
sualization  of  the  small  structures  for  the  sinogram 
without  averaging  while  the  visualization  is  about  the 
same  for  images  reconstructed  by  TV  and  FDK  from 
the  50-view  averaged  sinogram.  This  study  demon¬ 
strate  that  the  breast  C'T  dose  can  potentially  be  re¬ 
duced  to  a  very  low  level  for  screening  purpose  by  use 
of  very  sparse  projection  data. 


2 


